Multiphase coexistence in polydisperse colloidal mixtures 
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We study the phase behavior of mixtures of monodisperse colloidal spheres with a depletion agent 
which can have arbitrary shape and can possess a polydisperse size or shape distribution. In the low 
concentration limit, considered here, we can employ the free- volume theory and take the geometry of 
particles of the depletion agent into account within the framework of fundamental measure theory. 
We apply our approach to study the phase diagram of a mixture of (monodisperse) colloidal spheres 
and two polydisperse polymer components. By fine tuning the distribution of the polymer it is 
possible to construct a complex phase diagram which exhibits two stable critical points. 



I. INTRODUCTION 

In the study of the phase behavior of colloidal mix- 
tures there are two very distinct approaches. From a 
theoretical point of view one would like to keep the de- 
scription of the mixture as simple as possible in order 
to be able to calculate the phase diagram using well 
established theoretical tools. Furthermore, it is pos- 
sible to understand the main physics, i.e. the phase 
behavior, bulk correlation functions, or inhomogeneous 
density distributions, in simple terms. One example of 
such an theoretical approach is the Asakura-Oosawa-Vrij 
(AOV) model for colloid-polymer mixtures Within 
the AOV model, colloids are treated as hard spheres, 
while polymer behave like hard spheres only when inter- 
acting with colloids. The polymer-polymer interaction 
vanishes. Based on this model, one finds that the phase 
diagram of a colloid-polymer mixture, assuming perfectly 
monodisperse colloids and polymer, resembles that of a 
simple fluid, and the reservoir polymer density plays the 
role of an inverse temperature. One, relatively simple but 
successful, theoretical approach for treating the phase be- 
havior of this model colloid-polymer mixture is the free- 
volume theory (FVT) that assumes a low concentration 
of the depletion agent, the polymer. 4 

A very different picture of the phase behavior of 
colloidal mixture originates from experimental studies. 
While in some cases it is possible to prepare model sys- 
tems with almost monodisperse colloids and polymer f 5 -£ 
the generic case of a colloidal sample is polydisperse due 
to imperfections of the synthesis process. Polydispersity 
can be found both in size and in shape distributions. 
From the theoretical point of view, polydispersity is of- 
ten considered as a defect, because it complicates calcu- 
lations. However, polydispersity can give rise to inter- 
esting behavior, which can be utilized in experiments^ if 
the role of polydispersity is sufficiently well understood. 
One example is the study of the glass transition, where 
polydispersity in model hard colloidal spheres has the 
role as to suppress the freezing transition at high density 



of colloids so that a glass can form^ 

There are various very different ways to achieve a 
successful description. For hard and permeable spheres 
Salacuse and Stell 9 showed a thermodynamic approach 
based on the Percus-Yevick equation of state and the 
more accurate approach by Mansoori et al. The poly- 
dispersity in the chain-length of self-excluding polymer 
has been treated in Ref. [jjj- Polydispersity has been 
incorporated in the theory in different manner. Near- 
monodisperse distributions have been treated in a per- 
turbative fashionii and bimodal size distributions has 
been studied in Ref. [H. 

In the present study we wish to increase the under- 
standing of effects that can be caused by polydisper- 
sity further. To this end we consider a mixture of 
monodisperse colloidal spheres with a polydisperse de- 
pletion agent. We keep the theory simple by assuming 
that the depletion agent is efficient and can produce a 
sufficiently strong depletion attraction between colloidal 
spheres at rather low concentrations. In this limit we 
can generalize the free- volume theory within the frame- 
work of fundamental measure theory 13 (FMT) to cases 
of polydisperse size or shape distributions. 

As an application of our general and powerful approach 
we study the phase behavior of a ternary mixture of 
(monodisperse) colloidal spheres with two polydisperse 
polymer components with distinct size difference. Our 
study shows that the distribution of the depletion agent 
can be employed to tune the phase behavior of colloidal 
mixtures in a wide range and it is possible to generate 
a third fluid phase. Beside the low density (in colloids) 
gas and a high density liquid, we find a low density liquid 
phase. This is similar to findings in systems in which two 
length scales are introduced in the inter-particle interac- 
tion potentials ! 14 ! 15 ' 16 

The paper is structured as follows. In Section |TT] we re- 
call the theoretical basics of free-volume theory within 
the framework of FMT. We extend the approach of 
Ref. [13 to treat discrete and continuous distributions of 
the depletion agent. In Section HTT1 we discuss a few dis- 
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tribution functions which are frequently used in studies 
of polydisperse systems and show that for infinitely thin 
needles or platelets various types of size distributions lead 
to explicit expressions for the effective free-volume frac- 
tion a„Q. In Section IIVI we show and discuss as an ap- 
plication of our theory the phase diagrams for the case 
of a ternary mixture of (monodisperse) colloidal spheres 
and two polydisperse polymer components. We discuss 
the results and conclude in Section fVl 



II. THEORY 

We study mixtures of colloidal spheres immersed in a 
solvent of depletion agent. The latter component can 
possess a polydisperse size or shape distribution. Here 
we are mainly interested in mixtures where the depletion 
potential between the spheres is strongly attractive al- 
ready at rather low concentrations of the depletion agent. 
In this case we can employ the description of Asakura, 
Oosawa 1,2 and Vrij£ for the mutual particle interactions, 
which keeps the model simple. The colloid-colloid and 
the colloid-depletion agent interactions are modeled as 
hard-core interactions, while the mutual interaction be- 
tween particles of the depletion agent vanish. It has been 
shown 1 - that the phase behavior of the Asakura-Oosawa- 
Vrij model is closely related to the free- volume theory, in 
which the excess free energy of the mixture is expanded 
up to first order in the density of the depletion agent. 
This linearization of the excess free energy is equivalent 
to a vanishing interaction among particles of the deple- 
tion agent. As FMT— provides a good framework for 
the free-volume theory we first recall its main ideas as 
presented in Refs. [H and for discrete mixtures and 
generalize it to continuous distributions of the depletion 
agent. 



A. Free-volume theory 

In order to describe the phase behavior of colloidal 
mixtures one has first to choose the thermodynamic en- 
semble. For our purpose a system with a fixed number of 
big spherical colloidal particles N c , which is coupled to 
a reservoir of a depletion agent at chemical potential p d 
or fugacity Zd, respectively, represents the best choice.— 
Hence, the system is treated in a semi-grand canonical 
ensemble. The Helmholtz free energy separates into two 
terms 



f3F(N c ,V,z d ) = f3F(N c ,V,z d = 0) 



dz 



, r d(3F(N c ,V,z' d ) 



(1) 



The first term is the Helmholtz free energy of pure col- 
loids. By Taylor expansion of the integrand of the second 
term in Eq. {1} up to first order in z d we obtain an expres- 
sion for the volume available to particles of the depletion 



agent in a system of N c colloids at low concentrations of 
the depletion agent, the so-called free volume denoted by 
aV. In this limit of small fugacity of the depletion agent, 
we can rewrite Eq. (fT]) into 



0F(N c , V, z d ) = PF {N C , V) - p r d (z d )aV, 



(2) 



where p r d is the reservoir density of the depletion agent. 
Within the AOV model the free-volume fraction a de- 
pends only on the geometry of the depletion agent, e.g. 
the size ratio q = a d /a c in the case of spherical poly- 
mer, and on the packing fraction of the colloidal spheres 
Vc = ira 3 c N c /(6V). 

For the calculation of the phase diagram of the colloidal 
mixture we require the chemical potential of the colloid 
component p c and the osmotic pressure p. In the bulk 
two phases I and II coexist if they are in chemical and 
mechanical equilibrium, i.e. if = /i* 1 and = p , 
respectively. The chemical potential p c and the pressure 
p follow from the Helmholtz free energy via 20 



Pc 



dF 



9N n 



= Mo 07c) 



V,z d 



f da 



6 



(3) 



and 



P - 



dF 
W 



N c ,z d 



= Po(Vc) + P r d {z d ) [ a - r] c — ) . (4) 



The chemical potential and the pressure of the pure refer- 
ence system, labeled po and po, are described accurately 
by the Carnahan-Starling expressions in the fluid phase. 
To assure thermodynamic consistency the free-volume 
fraction a must also be calculated from the equation of 
state of the reference system, i.e. the Carnahan-Starling 
equation of stated Besides coexisting fluid states, i.e. 
points on the binodal line in phase diagram, we calcu- 
late the spinodal line which separates metastable states 
from thermodynamically unstable ones. For thermody- 
namically unstable states the second derivative of the 
Helmholtz free energy w.r.t. the colloid packing fraction 
•q c is negative and for points on the spinodal line it van- 
ishes, i.e. d 2 F/drjl = 0. 

For the solid-fluid coexistence we need the chemical 
potential and equation of state of a reference (fee) hard- 
sphere crystal, which is described by2i 



(3p = 2.1306 + 3 In 



3?/, 



Cp 



1]cp T]c / Vep Vc 



(•5) 



with rj cp = iry/2/6 ~ 0.74, the packing fraction for close 
packing, and^ 2 . 



f3pv s = 



3^c^?cp 
7]cp T]c 



(6) 



B. Fundamental measure theory 

For discrete mixtures of two species the excess 
Helmholtz free energy density <& = f3F/V is constructed 
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from two sets of weighted densities {n?} and {nf}, one 
for the colloids c and one for the depletion agent 
In homogeneous bulk fluids i labels four scalar-weighted 
densities, which can be identified as the scaled-particle 
theory variables Q of the iV-component mixtur o 23 ' 24 



N 

E^W 



(7) 



with the geometrical measures of component v. volume 
g v = v v , surface area gif 1 = a v , integrated mean curva- 
ture gv = c„ and Euler characteristics g^ = X v . 

An explicit expression for the free energy density of the 
AOV mixture in the limit of low density of the depletion 
agent can be obtained by expanding the hard-sphere (HS) 
mixture free energy with weighted densities n,i = nl + nf 
up to linear order in the density pd*^ 1 



AOV 



({<},{<}) = $ HS ({<}) 



3 

E 

k=0 



dni 



In order to relate the terms in Eq. ([5]) to the free-volume 
fraction a we transform the excess Helmholtz free en- 
ergy density, which in Eq. {T]) is given in the semi-grand 
canonical ensemble, into the canonical ensemble and ob- 
tain in terms of the colloid packing fraction r\ c and the 
density p d 



^{VcPdiq) = $0? c ) - Prfln(a(r? c ;g)). 



(9) 



By direct comparison of Eq. ([5]) and Eq. §§§ we obtain a 
FMT expression for the free-volume fraction 17 



a(Vc;{gl %) }) = exp 




gjKK}) W 
dni 9d 



(10) 



The partial derivatives of the free energy w.r.t. the 
weighted densities are thermodynamic quantities of the 
(pure) colloids, namely, the pressure po, the surface ten- 
sion 7 at a planar hard wall and the bending rigidities k 
and k, which describe the effect on the free energy due 
to curved surfaces i 17 i 25 ' 26 The depletion agent enters the 
expression for the free- volume fraction a only through its 

(k) 

geometrical measures g d . For the free energy density we 
employ Rosenfeld's original formulation of FMT 13 and 
the more accurate White Bear version from Rcfs. 27,28 



$({n a }) = -n ln(l - n 3 ) + 



1 - n 3 



+ 



n\ (n 3 + (l-n 3 ) 2 ln(l-n 3 )) 



367rn|(l 



(11) 



which is based on the Boublik-Mansoori-Carnahan- 
Starling-Leland (BMCSL) equation of state 2 ^ 



fipBMCSL = 



no 



71l71 2 



+ 



1-713 (1 - ti 3 ) 2 12tt(1 - 71 3 ) 3 

71 3 71 2 



36tt(1 - 7i 3 ) 3 



(12) 



It would also be possible to base the free-volume theory 
on a recent, slightly more consistent, mixture equation of 
state^" 3 ^ 



C. Generalization to polydisperse depletion agent 

FMT is a mixture theory by construction. For mix- 
tures of colloids and ideal depletion agent we can em- 
ploy parameters qi to parameterize the geometry of par- 
ticles of the depletion agent components. Although it is 
straightforward within FMT to treat general mixtures of 
v c colloid and Vd depletion agent components, we will in 
the following restrict our studies to the simpler case of a 
single colloid species with packing fraction rj c . For the de- 
pletion agent we start by considering i/j components. In 
the semi-grand canonical ensemble the excess Helmholtz 
free energy density is, analogous to Eq. ([9]), given by 



*(»7c, M<}; {*}) = * c 



Vd 

E 



4 a fe;{*})' (1 3 ) 



and the set {gi^} of geometrical measures is character- 
ized by the set of parameters {qi}. To study the influence 
of polydispersity in the depletion agent on the phase be- 
havior it might be sufficient to consider a discrete mix- 
ture, however we prefer to introduce continuous distri- 
butions, which can represent experimental systems more 
accurately than a discrete distribution. In the following 
we restrict the set {g„ } to depend on a single parameter 
q only. We introduce d(q) as a continuous distribution 
of the depletion agent in the reservoir, and require that 
J d(q)dq — 1. The density distribution of the depletion 
agent follows directly: p r d {q) = p r d d(q). By specifying the 
fundamental geometrical measures of the depletion agent 
as function of q we obtain, similar to Ref. 132 . 



Vc,p r d ; [d]) =®{Vc) 
-Pd, J d Q d (<l) a (Vc,Vd(q),ad(q),c d (q),X d (q)),(U) 

which is a functional of the distribution d(q). Note that 
the precise meaning of q is unspecified so far and can refer 
to a size ratio, in the case of spherical mixtures, or to a 
parameter that specifies the shape of the depletion agent 
in a more complicated manner. We return to this point 
later. Analogous to Eq. © we call the integral on the 
r.h.s. of Eq. (fl4| the effective free-volume fraction a e g 
and rewrite Eq. (fT4|) as 



PS 



$( Vc ) - p d a eS (r, c ;{d}). (15) 



III. DISTRIBUTION FUNCTIONS 

To study the phase behavior of a mixture of spherical 
colloids and a polydisperse depletion agent we focus on 
three frequently used distributions, namely the Schulz 
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(S), Gaussian (G) and the Hat-like (H) distributions. 
These distributions are characterized by two parameters 
corresponding to its first and second moment. The aver- 
age asymmetry ratio is given by q — &d/ ' <J Cl with the col- 
loid diameter a c and a length-scale of the depletion agent 
&d- The parameter z describes the degree of polydisper- 
sity: the limiting case z — > oo equals the monodisperse 
fluid and by decreasing the value of z the distribution is 
broadened. The distributions under consideration here 
lead to explicit expressions for the effective free-volume 
fraction a e g-(f? c ; d(q; g, z)), if the geometries of the deple- 
tion agent is simple. For more complicated shapes we 
have to perform the integration numerically. In terms of 
q and z the distributions considered here are given by 

MM,*) = (l)V- 1 ^^ , z> 1,(16) 
da{q;q,z) = -pexp(-(g- qfz 2 ) 

x2e(q)/[erf(zq) + l], (17) 

d H (q;q,z) = -^(q-q + z- 1 ) 

x6 (-q + q + z- 1 ) , qz>\. (18) 

It is evident that q > is required in all distributions and 
the integration in a e g-(ry c ; d(q; ; q, z)) is performed in the 
limits from to oo. In the case of a Gaussian distribution 
one can obtain simpler explicit expressions by imposing 
the full integration range q = — oo to oo. The result, d' G , 
is similar to that of the cut Gaussian distribution da 

d' G («?;g,z) = ^=exp(-(g-g) 2 z 2 ). (19) 

Obviously, it is possible to consider different size distri- 
butions such as the log-normal distribution, which decays 
slower for large values of q than the distributions consid- 
ered here. 



A. Explicit expressions for Q e ff 

The equations for the chemical potential /x c and the 
osmotic pressure p from Eqs. ([3]) and (|4|) for polydis- 
perse size distributions contain integrals over the distri- 
bution. For the calculation of the phase diagram these 
equations must be solved simultaneously. In general 
this must be done numerically In the case of one- and 
two-dimensional depletion agents, such as infinitely thin 
rods or platelets, their geometry is sufficiently simple so 
that the effective free-volume fraction a e g can be calcu- 
lated explicitly. Infinitely thin needles represent an one- 
dimensional depletion agent. Infinitely thin platelets rep- 
resent a two-dimensional depletion agent. Here we con- 
sider both disk- and hexagon-shaped platelets. The geo- 
metric measures for these geometries are given in Tab. HI 
In Appendix we show explicit expressions for a e ^ for 
one- and two-dimensional depletion agents and the dis- 
tributions given above. As aforementioned, we obtain 



Dim Depletion a; 


;ent 


Geometry 


I'd 


CLd 


Cd 


x d 


1 needle 


L — q a c , R = 








L 
4 


1 


2 disk 




R = qa c /2, h = 





2tvR z 


\R 


1 


2 hexagon 




R = qa c /2, h = 





Z-JlR 2 




1 



TABLE I: Geometrical measures for one- and two-dimensional 
depletion agents. All lengths are measured in units of the 
diameter <r c of the spherical colloids. We consider needles 
of length L and radius R = and platelets with radius R 
and thickness h — 0. The value of q plays the role of a size 
ratio and compares the size of the depletion agent to that of 
a colloid. 

the thermodynamic quantities of the pure hard-sphere 
fluid that enter the free-volume fraction a e g from the 
excess free energy density $ of the White Bear version 
of These quantities are the pressure (3p — 

<9$/<9n 3 , the surface tension /?7 = d^/dn 2 and bending 
rigidities /3k — <9$/<9ni and j3R = d^/dn ^^^ 

In the present study we restrict the product qz to suf- 
ficiently large values. For the Hat-like distribution we 
require qz > 1 which ensures that J °° dqd(q) = 1. In the 
case of the Gaussian distribution d' G we require qz > 1 
because for small values of qz the contribution of the dis- 
tribution from values q<0 increases, e.g. for q — 0.5 and 
z = 3 we obtain J^dqddq) - 0(1CT 2 ). 

In the case of two-dimensional depletion agents, the 
geometric measures can differ for different geometry of 
the depletion agent: for disk- and hexagonal platelets ad 
and Cd differ - see Tab. HI For that reason we show the 
result for a g g: in Appendix [X] in its general form. 

IV. RESULTS 

While our approach developed in Section |TT] is gen- 
eral and allows one to study various systems, we want 
to present only one application. We study the phase 
behavior of a ternary mixture of (monodisperse) col- 
loidal spheres with two polydisperse polymer compo- 
nents of distinct sizes. Both polymer components pos- 
sess a size distribution which we describe by Schulz dis- 
tributions, which in total lead to ds(q;q,Q, Zg, zq) = 
x ds(q; q, Zg) + (1 — x) d$(q; Q, zq) with a mixing parame- 
ter x and two averaged size ratios q and Q. For simplicity 
we assume the width of both parts of the distribution to 
be equal. The polymer reservoir packing fraction is given 

by v r d = f (q) 3cr cPd> where (i) = x $ + (i - x ) Q- 

We expect the phase behavior to be most interesting if 
the average size ratios q and Q are distinct, so that the 
phase diagrams in the limiting cases x = and x = 1 
are sufficiently different. To this end we choose q = 0.25 
and Q = 2.0. In the case of x = when we obtain a 
binary mixture with an averaged size ratio Q = 2.0, the 
phase diagram exhibits a stable fluid-fluid phase sepa- 
ration into a colloid poor gas and a colloid rich liquid 
phase. For the mixing parameter x = 1, when the av- 
eraged size ratio q = 0.25, the fluid-fluid phase separa- 
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tion, which is metastable w.r.t. to crystallization in the 
monodisperse case^ can be stabilized by a sufficient de- 
gree of poly disp er sity . 33 > 34 ' 3 5 

Before we study the full phase diagram of the ternary 
mixture, we determine which value of Zq is required in 
order to stabilize the fluid-fluid phase separation for x — 
1. 



A. Influence of size polydispersity 



n T3 
Q_ 
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In this section we study the phase diagram of a binary 
colloid-polymer mixture, which corresponds to the limit- 
ing case x = 1, with an averaged size ratio of q = 0.25. 
In contrast to the aforementioned examples of one- and 
two-dimensional depletion agents, a e g can not be ob- 
tained explicitly. Therefore, we have to determine the 
phase coexistence by solving two equations, Eqs. ([3]) and 
(|4]), numerically. 

The equation of state of the pure reference system and 
the free-volume fraction a e fj respectively, and hence the 
chemical potential fi c and the osmotic pressure p of the 
colloids, are based on the free energy density <I> of the 
White Bear version of FMT, Eq. (fTT|) . 

In Fig. QJa) we show the phase diagrams of polydis- 
perse AOV mixtures as we vary the degree of polydisper- 
sity. For the size ratio q = 0.25 and a narrow distribution 
of the polymer, the fluid-fluid phase separation remains 
metastable similar to the monodisperse case. We show 
the phase diagram for z = 50 (full lines) in Fig. [lja) and 
the corresponding size distribution in Fig.[ljb). Upon de- 
creasing z, which broadens the size distribution, the fluid- 
fluid and the fluid-solid coexistence lines come closer to 
each other and almost touch for z — 10 (not shown). By 
further increasing the degree of polydispersity the fluid- 
fluid coexistence becomes stable, as we show in Fig. Ufa) 
for z = 5 (dashed-dotted lines) and z — 2 (dotted 
lines). The corresponding size distributions are shown 
in Fig. |ljb) . The critical point is shifted to lower colloid 
packing fractions rj c , as the value of z decreases. 

The polymer distribution in the system differs from 
that in the reservoir^ which is a phenomenon known 
as fractionation. For spherical polymer we find results 
equivalent to those of Ref. |3J. At low packing fractions 
of the colloids the size distribution of (spherical) polymer 
follows closely the reservoir distribution. As the colloid 
packing fraction increases, the maximum of the poly- 
mer size distribution decreases significantly and moves 
towards smaller polymer radii. In principle we observe 
a qualitatively similar but less pronounced behavior also 
in the case of platelet-like or rod-like depletion agents. 
The maximum in the size distribution of the depletion 
agent becomes smaller and moves to smaller sizes as the 
colloid packing fraction in the system increases. This is 
to be expected, because one- or two-dimensional deple- 
tion agents require less free volume to fit between colloids 
than spherical polymer. 

We have confirmed that a treatment of the phase be- 
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FIG. 1: (a) Phase diagrams of mixtures of colloids and 
polydisperse (spherical) polymer with an average size ratio 
q — 0.25 and (b) the corresponding distributions. We em- 
ploy a free-volume theory based on the BMCSL equation 
of state. For low degrees of polydispersity, corresponding 
to z — 50 (full lines), the fluid-fluid phase separation is 
metastable w.r.t. fluid-solid coexistence. Upon increasing the 
degree of polydispersity, corresponding to z = 5 (dashed- 
dotted lines) and z = 2 (dotted lines), the fluid-fluid coex- 
istence is stabilized. 



havior completely based on the PY compressibility equa- 
tion of state yield qualitatively similar results. In the 
present approach this is achieved by employing the free 
energy density $ of Rosenfeld's original formulation of 
FMT to evaluate a e g and the pressure and chemical 
potential for the pure reference system. The results in 
Refs.HHl and k35l are somewhat in between these two 
predictions as their approach employs a hard-sphere ref- 
erence system based on the Carnahan-Starling equation 
of state and a free- volume fraction based on the PY com- 
pressibility equation of state, which for spheres is equiv- 
alent to the scaled-particle theory expression. 
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B. Influence of shape polydispersity 

In this section we briefly address the question whether 
polydispersity, which describes a shape distribution, can 
also stabilize the fluid-fluid phase separation. One can 
think of the shape polydispersity in two different ways. 
On the one hand it can be thought as describing effects 
due to imperfections in the synthesis of colloidal parti- 
cles. On the other hand, and this is more relevant in the 
present context, it can be seen as a simple model for the 
effect that polymer coils are not rigid and can change 
their average shape when squeezed in between colloids. 

To this end we consider ellipsoids, with half-axes 
a = b = <Jd/(2y/q) and c = qad/2, for which the ge- 
ometrical measures can be obtained explicitly^ This 
case is relatively simple as a single parameter q char- 
acterizes the shape while ad — 0.25cr c is fixed. Ob- 
viously, more complicated scenarios are possible. Note 
that for this choice of half- axes the volume Vd = 
is kept fixed, while the surface area = ira~i(l/q + 
q 3 / 2 (arccosh l/q 3 ^ 2 )/\/l/q — q 2 )/2 and integrated mean 
curvature Cd — &d(q + (arccosg 3 / 2 )/ \J q — q 4 )/<i depend 
on q. The Euler characteristic Xd — 1, independent of 
q. Here, the parameter q describes the degree of devi- 
ation of the shape from a sphere: for small values of q 
we obtain lens-shaped particles (oblates), while for large 
values of q the particles become cigar-shaped (prolates). 
Similarly to the case of size polydispersity, we observe 
that an increasing degree of shape polydispersity stabi- 
lizes the fluid-fluid phase separation w.r.t. fluid-solid de- 
composition. In Fig. [2] we show an example for a narrow 
shape distribution with z = 50 (full line), in which case 
the fluid-fluid critical point is metastable, and one for a 
broad distribution with z = 2 (dotted line). For z — 2 
we find a stable fluid-fluid phase separation. 

It seems that polydispersity in general favors a stable 
fluid-fluid phase separation because the polydisperse de- 
pletion agent can fill the free volume of the system more 
effectively than a monodisperse one. Even if the bigger 
particles of the distribution cannot find free volume in the 
system, the smaller ones still do. In terms of the deple- 
tion potential this would result in a longer range of the ef- 
fective attraction between colloids ^ giving rise to a more 
negative effective second virial coefficient Bf which in- 
creases the tendency for fluid-fluid phase separation; 37 
This effect seems to be robust against details of intro- 
ducing polydispersity and can also be observed for a de- 
pletion agent with a = b = erd/2 and c = qad/2^ where all 
geometrical measures, except of the Euler characteristics, 
depend on q. 



C. Ternary mixtures with a bimodal polydisperse 
distribution 

Based on our findings in Section TlV A I we choose Zq — 
zq = 5 for the bimodal polydisperse polymer distribu- 
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FIG. 2: Effect of shape polydispersity of the depletion agent 
on the phase behavior of a mixture of colloids and depletion 
agent. The depletion agent consists of ellipsoids with half- 
axes a = b = ad/(2y/q) and c = qad/2. The mean value of 
the distribution is fixed at q — 1.0 while the size ratio between 
the depletion agent and the colloid is fixed at ad/cr c = 0.25. 
For 2 = 50 (full lines) the fluid-fluid coexistence is metastable 
and similar to the monodisperse case (not shown here). Upon 
increasing the degree of polydispersity we find two stable fluid 
phases. The phase diagram for z — 2 (dashed lines) exhibits 
a stable critical point and a triple point. 

tion. 

Before we show the calculated phase diagrams, it is in- 
teresting to compare the two stable critical points corre- 
sponding to the limiting cases, i.e. x = 1 and x = 0. The 
critical points can be found at very distinct regions of the 
phase diagram. In Fig. [3] we show the paths of the criti- 
cal points through the phase diagram as a function of the 
mixing parameter x. When we start at x = 1, for which 
only the smaller polymer are present, we find a critical 
point at rj c = 0.258 and rf d = 0.247. As we decrease 
the value of x, the critical point moves towards smaller 
values of rj c until this critical point becomes metastable. 
For x < 0.979 the critical point no longer exists. When 
we start at x — 0, for which only the larger polymer are 
present, we find a second critical point at r\ c = 0.014 and 
rf d = 0.924. As the value of x is increased we find that 
this critical point moves towards larger values of rj c until 
it becomes metastable. For x > 0.997 this critical point 
vanishes. 

It is most interesting to note that both critical 
points found in this system, which correspond to criti- 
cal points of the limiting cases, can exist at the same 
time for a range of values of the mixing parameter, 
x e [0.979,0.997]. In order to obtain a full picture and 
determine where these critical points are stable we cal- 
culate the full phase diagrams in this regime. 

In Figs. [UJa)-(b) we show phase diagrams in the r\ c - 
r72 - representation for different values of x. For x = 
0.9999 [see Fig. Ufa)] the phase diagram is similar to the 
monodisperse case studied above. Upon decreasing x, the 
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fraction of polymer with average size ratio Q increases, 
and the spinodal displays two minima, one correspond- 
ing to the stable critical point (triangle) and the sec- 
ond minimum corresponding to an additional metastable 
fluid- fluid phase separation (diamond) - see Fig. HJb). 
Further decreasing the value of the mixing parameter to 
x = 0.9935 [Fig. HJc)] yield a phase diagram with two 
stable critical points and a second triple point. At one 
triple point, a gas, a liquid, and a solid phase coexist, 
while at the second triple point a gas phase and two liquid 
phases with different densities coexist. There is a novel 
low density liquid phase in the region between the two 
critical points with colloid packing fractions r\ c = 0.0177 
and 0.251. As we vary the value of x the triple points 
move. According to Gibbs' phase rule a four-phase co- 
existence between a gas, the low density liquid, the high 
density liquid and the solid is also possible. For our sys- 
tem it would be expected for x between the values used 
in Fig. [Kb) and (c). In Fig. GJd) for x = 0.98 we also 
observe two critical points. The critical point with lower 
colloid packing fraction (diamond) is stable. The second 
critical point at higher colloid packing fraction (triangle) 
is now in the metastable regime and vanishes upon de- 
creasing the values of x. For these systems one liquid 
phase and one triple point are stable. Further decreasing 
the values of x below approximately 0.8 lead to a shift of 
the coexistence region to higher values of ff d . For x <C 1 
the fluid-fluid binodal is in the region of the fluid with 
monomodal polydisperse depletion agent with Q = 2.0. 

We observe similar behavior for a mixture of colloids 
and two monodisperse polymer components, i.e. z q and 
zq — > oo. For the same size asymmetries as above, 
q — > q = 0.25 and Q — > Q — 2.0, the critical point is 
metastable when x — 1, or equivalently (q) = 0.25. We 
observe one stable and one metastable critical point for 
low colloid packing fraction upon decreasing the value of 
x. In the case q = 0.5 and Q = 2.0 there is a stable 
critical point for x = 1. Decreasing x, similarly to the 
polydisperse case described above, leads to a second sta- 
ble critical point. To summarize, we find that both by 
means of choosing the degree of polydispersity and the 
size asymmetry we can treat a model fluid which exhibits 
one or two (stable) critical points and one or two (sta- 
ble) liquid phases. Furthermore, we have verified that it 
is possible to generate two critical points and hence three 
fluid phases by using a depletion agent of different than 
spherical shape. 



V. DISCUSSION 

We have generalized a novel approach to free-volume 
theoryil based on FM T 13 i 23 ' 24 to treat a depletion agent 
of arbitrarily shaped particles with a continuous distri- 
bution. Here, we introduced one parameter, q, to specify 
a continuously varying distribution of depletion agent. 

To evaluate the phase diagram within FVT an expres- 
sion for the free- volume fraction is required. In the case of 
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FIG. 3: The path of the two critical points in a ternary mix- 
ture of colloidal spheres and two polydisperse polymer com- 
ponents. If the mixing parameter x = 1, only the smaller 
polymer are present and we find a stable critical point at 
r/ c — 0.258 and r/2 = 0.247. As we decrease the value of x this 
critical point moves towards smaller values of r\ c and vanishes 
at x « 0.979 (upper line). For a mixing parameter x = 0, 
when only the bigger polymer are present, we find a second 
critical point at r\ c = 0.014 and r[ p = 0.924, which moves to- 
wards larger values of 77c (lower line). For x G [0.979,0.997] 
(indicated by the dotted lines) we find two critical points in 
the system. 



polydisperse size or shape distributions for the depletion 
agent, one can obtain an effective free-volume fraction, 
which in general has to be calculated numerically. In the 
cases of infinitely thin needles or platelets it is possible to 
get explicit expressions for the effective free- volume frac- 
tion for various distributions which are of interest from 
an experimental and a theoretical point of view. We give 
some expressions for a e g- in Appendix [Al 

Beside size polydispersity, we have also considered 
shape polydispersity. To this end we have calculated 
phase diagrams of a mixture of colloidal spheres with 
penetrable ellipsoids with half-axes a = b = (Td/(2^/q) 
and c = qtJd/2- 

Our results confirm that polydispersity has the ten- 
dency to stabilize the fluid-fluid coexistence region 
w.r.t. crystallization. We have discussed this trend, in 
agreement with earlier studies , 33 i 34 i 35 in Section HV Al for 
size polydisperse spherical polymer coils. In Section flVBI 
we have shown that similar behavior can be found for 
the case of shape polydisperse case. Upon increasing the 
degree of polydispersity the fluid-fluid binodal "crosses" 
the fluid-solid coexistence line and one obtains a stable 
critical point with associated liquid and gas phases and 
a triple point. This effect seems to be robust against 
various ways to incorporate polydispersity, which can be 
rationalized by the fact that particles of a polydisperse 
depletion agent can fill the available free volume in a sys- 
tem better than a monodisperse species. 

Our main result is the study of mixtures of colloids and 
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FIG. 4: Phase diagrams for a mixture of colloids and a bimodal distribution of polydisperse (spherical) polymer. For x very 
close to 1 we recover the phase diagram for the monomodal distribution, see plot (a) for x = 0.9999. The fluid-fluid coexistence 
(full line) is stabilized by polydispersity. Upon decreasing x, i.e. by adding larger polymer, the spinodal (dotted line) exhibits 
a second minimum for x = 0.995 in Fig. (b) and we observe a second critical point, which is metastable (diamond). For 
£=0.9935 we observe two stable critical points. Beside a low density (in colloids) gas (G) there exists a low density (LDL) and 
a high density (HDL) liquid phase - see Fig. (c). For a value of x — 0.98, shown in (d), the critical point at high values of rj c 
(triangle) is metastable. 



two species of polydisperse (spherical) polymer coils. The 
polymer were described by a bimodal distribution con- 
sisting of two polydisperse distributions. For the aver- 
aged size ratios of the polymer components chosen here, 
i.e. q = 0.25 and Q = 2.0, we find in distinct parts 
of the phase diagram fluid-fluid phase separations in the 
two limiting cases x = 1, when only the smaller poly- 
mer species is present, and x — 0, when only the larger 
polymer species is present. For x = the correspond- 
ing fluid-fluid phase separation is stable independent of 
the degree of polydispersity. In the second limiting case, 
x = 1, the fluid-fluid phase separation has to be stabi- 
lized by a sufficient degree of polydispersity. We choose 
Zq — zq — 5 and find for a range of the mixing param- 
eter x two stable liquid phases: beside the gas phase, 
which is diluted in the density of colloids, we find a low 
density and a high density liquid phase. The fluid-fluid 
coexistence region is w-shaped in this case and the phase 



diagram exhibits two stable critical points and two triple 
points. 

For systems for which the inter-particle interactions 
contains two length scales, such as short-range repulsion 
and long-range attraction, Hemmer and Stell 14 ' 15 found 
that two stable critical points may occur. More directly 
related to our situation is a system of star polymers 16 
in which similar features in the phase diagrams are re- 
ported, as the effective interaction potential is compa- 
rable to those of the model of Stell and Hemmer. For 
other model fluids similar effects were observed i 38 ' 39 ! 40 
We point out that our findings are not restricted to a 
polydisperse depletion agent. We find similar results for 
a ternary mixture in which we consider one colloid species 
and two distinct monodisperse species of depletion agent. 

While often size or shape polydispersity is consid- 
ered an undesirable side-effect, originating from imper- 
fections of particle synthesis, it is possible to make use 
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of this property by stabilizing fluid-fluid phase separa- 
tions. Moreover, the degree of polydispersity allows one 
to adjust the phase diagram. From this point of view 
one might speculate about fine tuning the phase behav- 
ior of colloidal mixtures by specifying the polydisperse 
distribution of the depletion agent. 



APPENDIX A: EXPLICIT EXPRESSIONS FOR 

Q eff 

Explicit expressions for the effective free-volume frac- 
tion in the case of infinitely thin needles for the Schulz, 
Gaussian and Hat-like distributions, Eqs. (fl6|) - (fl"8l) . read 



a e Q G (r] c ;q,z) = exp I —k - qnR c /2 + 



a eS,s(^*) = {9KR *+ 2z )'e* > (A1) 

16z 2 

/ k 2 R 2 \ 

a 'eS,G^ c '^' z ^ = exp ( ~ R ~ V kRc / 2 + )' ( A3 ) 

2z /kR \ 

a e^_,H^^ z ) = ^R~ exp(_S ~ qnR c /2)smh ( j . (A4) 

I 

Since only k and R enter these expressions for a one- tion of state. In terms of the colloid packing fraction rjc, 
dimensional depletion agent, there is no difference be- average asymmetry ratio q and the polydispersity param- 
tween an approach based on the PY or the BMCSL equa- eter z we obtain 



a eSs (r) c ;q,z) = (2z)*(l-r? c ) (^^- + 2zj , (A5) 

' ^l-erff-zg+ - X ], (A6) 



eii(zq) + 1 V V 4z (! - Vc) 

2(1 -Vc) 16z 2 (l-, ]c ) 



a'esJVc; q, z) = (1 - ry c ) exp - _ c + e , (A7) 



_ 2z(l- Vc ) 2 ( -3q Vc \ . ( 3r; c \ 

9, *) = 3??c ) «*h ^- (1 - ^ j . (A8) 

I 

In the case of infinitely thin disk- or hexagon-shaped fraction within the BMCSL approach are given by 
platelets the expressions for the effective free-volume 



OeEsfocjg,*) = (I) (-2 X2 r z/2 exp(-|i+ Xo )Z}_ 



8X2 



-Xi 



(A9) 



a eff g(^'9i z ) = ex P 



z 2 q 2 X2 + Xo z 2 - Xo X2 + z 2 qxi + Xi 2 /4 
z 2 - X2 
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x erf 



z 2 g + xi/2 



a 'eS,G^ c;qiZ ^ = exp 



\/z 2 - X2 J J 2[erf(zq) + l] v /z 2 -X2 
z 2 q 2 X2 +Xoz 2 -X0X2 + z 2 qxi + Xi 2 / 4 \ z 



Z - X2 



- X2 ' 



«eff,H(^;3^) = i^/^ ex p(xo-^) ^(-l) n erf 



xMyv iNn„„ r / /2 X2(9-(-i) n ) + xiz 

n=l 

I — 



(A10) 
(All) 
(A12) 



where X2 = —PlQ-d/q 2 , Xi = —Pucd/q and xo = —finXd- Tab.|TJ D_ z with z > are parabolic cylinder functions. 
The geometrical measures a^, and A^ are given in 
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